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Abstract 

We discuss how magnetic screening can be systematically included in a self-consistent 
way for Chromodynamics at high temperatures. The resulting gap equation, which gives 
the summation of self-energy insertions, is calculated to one-loop order and leads to a 
nonzero value for the magnetic mass. 
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The generating functional for hard thermal loops in Quantum Chromodynamics (QCD) 
is closely related to the eikonal for a Chern-Simons theory [1] . The importance and various 
properties of hard thermal loops have been the subject of many recent investigations [2-5] . 
The generalization of the Chern-Simons eikonal to a moving plasma suggests that there is 
another closely related gauge-invariant mass term which gives screening of static magnetic 
interactions, in other words, a magnetic mass term [6]. Although nonlocal, the equations 
of motion for this term can be written in a local way by using auxiliary fields and without 
introducing additional degrees of freedom. Because of this last property, the ultraviolet 
behaviour of the theory is unchanged. On this basis, it was suggested that this mass term 
could be used as a gauge-invariant infrared cutoff for loop calculations in QCD at high 
temperatures. In this paper we discuss how this can be carried out systematically. Some 
recent related papers are listed as ref. [7]. 

A consistent implementation of a gauge-invariant infrared cutoff will lead to a gap 
equation for the value of the magnetic mass. We obtain this equation to one-loop order; 
as usual, this means a self-consistent summation of one loop self-energy insertions. The 
magnetic mass to this order is obtained as ^.ZdiCg^T / Att where g is coupling constant, T 
is the temperature and C is defined by C5°'^ = fo-mnj^bmn^ jamn ^ei^g ^j^g structure con- 
stants of the gauge group. (C = N for an SU{N)-gQMge theory.) An immediate question 
is whether the two- and higher loop contributions are smaller than the one- loop terms. 
Purely based on counting dimensions of momentum integrals and powers of , we cannot 
conclude whether higher loop effects are smaller or not. It is really an issue of numerical 
factors and possible logarithms of g. Now, the magnetic mass, by the standard arguments 
of dimensional reduction at high temperatures is related to the mass gap of the three- 
dimensional QCD. For the mass gap of three-dimensional gauge theories, the one-loop 
calculations can hardly be adequate. However, for the quark-gluon plasma, perturbative 
loop expansion, perhaps with resummations, is expected to be reasonable at high tem- 
peratures, at least for a number of processes of interest and our approach is applicable. 
In any case, qualitatively, it is interesting that a nonzero value can be obtained to this 
order; computationally, it is useful to have a loopwise implementation of gauge-invariant 
magnetic screening, irrespective of the specific numerical value. 

In the rest frame of the plasma, the magnetic mass term considered in [6] has the 
following form 

f = -M^Sm (la) 
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Sm^ / dnK{Ar,,Ar,) (lb) 



where = ^AiTii^An = ^Aifti. rii is a (complex) three-dimensional null vector of the 
form 

rii = (— cos 6 cos (p — ism(p, — cos 6sm(p + i cos (p, sin 6) (2) 

In Eq.(l), dQ — sin 6* d9d(p and denotes integration over the angles of nj. K{An,Afi) is 
given by 

K{Ar^,Afi) = ~ J d^X^ J dh Tr{Ar^,Ar^) + i7^I{Ar^) + i7^I{An) 

z = n ■ X, z = n ■ X, x'^ denotes coordinates transverse to n, i.e., x^ ■ n = 0. Also 



(3) 
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Zij = Zi — Zj. The argument of all A^s in Eq.(4) is the same for the transverse coordinates 
x^. The lowest order term in Sm was shown to be 

JwT' ^^^^-^^^^"^^^ {^'^ ^^^'^ 

This term involves only the transverse potentials, as expected for magnetic screening and 
on account of gauge invariance. The terms with higher powers of A make 5'^^ invariant 
under the full non-Abelian gauge transformations. 

The general strategy for the inclusion of this term is as follows. We write the action 

as 

S^So + M^Sm -ASm (6) 

^0 is the standard quark and gluon part of the action. Below we shall not consider the quark 
terms since their effects are small and have the same general qualitative features. Sq will 

be just the Yang-Mills action. A is taken to have a loop expansion, A = A^-*^) + A^^) + 

Calculations can be done in a loop expansion. We require the pole of the propagator to 
remain at kQ — = (for the transverse potentials) as loop corrections are added. 
This requires choosing A^-*^) to cancel the one-loop shift of the pole, A^^) to cancel the 
two-loop shift of the pole, etc., as is usually done for mass renormalization. Of course, 
we do not want to change the theory, only rearrange and resum various terms. Thus we 
should impose the condition 

A = A(i) + A(2) + . . . = (7) 
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This condition is the gap equation determining M in terms of g'^ and T. This procedure of 
adding and subtracting a mass term, with a gap equation required for consistency, is very- 
standard, for example, in the Nambu-Jona Lasinio model. It amounts to the self-consistent 
summation of self-energy corrections. The difference in the present case is that, for reasons 
of gauge invariance, the mass term involves an infinite number of interaction vertices as 
well. 

We shall now turn to the explicit one-loop calculations. For the plasma at high 
temperatures, the effects of magnetic screening will be insignificant for high momentum 
processes. The regime of interest involves momenta small compared to T. The thermal 
part of the gluon propagator simplifies as 

l^,/T _ 1 ^ ^ ['^(^0 - cufc) + 5(A;o + uju)] (8) 

ujk = -h M^. This is equivalent to using Euclidean three-dimensional propagators, 
with a coupling constant e — \J g'^T] the thermal part of the loop contributions can be 
done in a three-dimensional theory. This is, of course, the standard dimensional reduction 
argument. The electrostatic field, with a Debye mass of order ^T, will also be neglected for 
low momentum calculations. The relevant momenta for which this approximation is valid 
will be of order ^r^T; the coupling constant is also self-consistently evaluated at a scale of 
order g'^T. This is all in keeping with the assumed hierarchy of Aqcd « g^T « gT « 
T for the hot quark-gluon plasma. 

The action for momenta small compared to T and gT can be written as 

A 



1 

rpa rpa , n 



^ / d-'x (9) 



S = j d^x 

F^^ = d^Al - dyAl + f'^'^A^^Al. The integral of Cm is the thre^dimensional Euclidean 
version of Eqs.(lb,3,4). There is now only one transverse coordinate x'^ . A convenient 
gauge-fixing term is |5 • ^(1 — M'^^)d ■ A. This gives the gluon propagator as S'^^Sij {k^ + 
M'^)~^. The sum of all one-loop contributions to the gluon polarization is finite. The 
functional integral is thus 

J [dA] det{-d ■ D)e~^ (10) 
D^^ = 8^5"^ + Z'*'^''^^. The action simplifies as 



S = j |^«(-a2 + M2)^« + j A-A'jAir'^Vijk + Vijk) 

I 1_ jamn jabc j J^. A^. -|- 



(11) 



Vijk{k, q, -{k + g)) = - [{2k + q)jdik - {2q + k)i5jk + {q- fc)fc%] 



Vijk{k,q,-{k + q)) = 



M2 

24^ 



k ■ n 



q ■ n {q + k) ■ n 
q ■ n {q + k) ■ n 



(12a) 
(126) 



This, after angular integration, becomes 



Vijk{k,q,-{k + q)) = -i 



1 



k-{q + k) 



q ■ k q ■ (q + k) 



kikjkk 



+ 



{q + ky 



B {q + ky 
{qiQjkk + qkQikj + qjqkh) - {q ^ k) 



(12c) 

Expression (12c) is the contribution from Sm- The four-point vertex involves, in addition 
to the standard Yang-Mills vertex of Eq.(ll), a term from Sm of the form J dfl (ninjUkni). 
Since Hi is a null vector we get zero from this term to one loop order. This is one of the 
advantages of our form of Sm', because rii is null. Wick contractions at the same point, 
with a propagator of the form 6ij{k'^ + M^)~^, give zero. Vertices of higher than quartic 
order do not contribute at one-loop level. The relevant Feynman diagrams are shown in fig. 
1. The calculations are straightforward although somewhat tedious. Despite the algebraic 
complexity of Vijk, there are many cancellations and for the one loop part of the effective 
action we find 



^ '^At{-k)U,j{k)A^{k) + 



n 



(13) 
(14a) 



n 



(1) 

^3 



C5. 
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a — a^ 



+c 



c_ 
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Jvj 

IGtt 

1 



da 



167r ^k'^{a-a'^) + M^ 

3 + 6(a-a2) 2 /- ^' 

y/k'^{a-a^) + M^ k^^ \ 



(146) 







k"^ 



da 



da 



\/B{a - a^) + M2 - ^k'^{a - a^) + M^a 
2yjk'^{a - a^) + M2q; - ^Jk'^{a - a^) + M2 

—k\/ a — 



(14c) 
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n^j-* is the contribution of the standard Yang-Mills diagrams, n^^-* involves the new vertex 
(12c), by itself and mixed with (12a). The a- integrals can actually be evaluated in terms 
of elementary functions but we do not need the explicit form in what follows. The total 



contribution to Ily can be written as 



n 



kikj 



(15a) 



U{K) 



47tM 



da 



^^K{a-a'^) + l-a + a^ + 1 - 2y/K{a - a'^) + a'^ 



(K-l) 



(5 - 2a + 2a^) 



4 ^/K{a - a'^) + 1 - a + a'^ 



-bVl-a + a^ (15a) 



B 



da 



— a-\- + 



4 Vl-a + a2 



4 Vl - a + a2 



-3 



(156) 



: (2.384) 



47rM 

K = (A;^ + M'^)/M'^. Notice that Ily is transverse as expected on grounds of gauge 
invariance. 11 (K) vanishes at i^T = and has no other zeros for positive K as can be checked 
graphically. (For negative K, n(i^) is complex and our calculation which uses time-ordered 
products is not applicable; see ref. [4].) The one- loop corrected inverse propagator has 
the form ^[K + Ii{K) + {B - ^)] . Since U{K) = at K = 0, we see that the pole of 



the propagator will not be shifted if we choose B 
equation (7) then gives 



M2 



or 



4tv 



M ^ (2.384)-=— ^ (2.384)^ 



(2.384). The gap 



(16) 



47r Att 

With this choice of A^^\ the inverse propag ator has the form ^[K + U{K)] . The correc- 
tion n(i^) is significant compared to K, for small K; for large K, it approaches the result 
for massless Yang-Mills theory; this is a necessary check for any infrared cutoff. 

Eventhough Ily is transverse in agreement with gauge invariance, it depends on the 
gauge fixing used for the gluon propagator in the loop. The position of the pole and hence 
the gap equation do not depend on this. One can explicitly check this. The simplest way 
is as follows. In the effective action, one can have a gauge-dependent 11^; Generally, the 
higher point functions are gauge dependent as well. Physical results, such as scattering 
amplitudes, are independent of the gauge fixing used. Alternatively, one can define a new 
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two-point function, a new three-point function, etc., which are independent of the gauge- 
fixing and which lead to the same physical results. One shifts some of the contribution 
(to the scattering process) from the three-point vertex to the propagator; similar shifts are 
done for the higher point functions as well. The amount of shift is determined by requiring 
the physical results to be the same and leads to the pinching procedure [8] . For our gauge 
choice, the pinching terms arise from the diagrams shown in fig. 2. (A quark scattering 
process suffices to identify these terms.) 




'WWWV^ 

\ 





New vertex 



Standard YM 



figl.- Gluon self-energy 





fig2.- Pinching diagrams, 



The total pinching contribution to the two-point function is 

U{K) = f da 

^ ' 167r io 



2a 



, , (17) 

^/K{a^^oF)Tl^^^^To? y/Kia - a^) + 

This is to be added to n(i^) to obtain the gauge-independent two-point vertex function 



From the explicit factor of K in Eq.(17), we see that the pole remains at K = 0. 

We shall now briefly consider two-loop corrections. These must involve and since 
has the dimension of mass, the two-loop contribution to B has the form 

B{2-loop) ^ ^ log(T/M)] (18) 

where 71,72 are pure numbers. The gap equation now reads 

a;' = ^ X + 71 + 72 log(T/M) (19) 

M = Ce^x. The numerical values of 71,72 determine whether the two-loop corrections are 
small. There are seventy-nine diagrams for 11^ at two-loop level as opposed to five at the 
one-loop level. Since the vertices given by Sm are fairly complicated we have not completed 
the two-loop calculation. A preliminary analysis of some of the diagrams suggests that 
the two-loop effects may be smaller by a factor of three or four. It should be noted that 
this issue does not impinge on the use of Sm as a gauge-invariant infrared cutoff using the 
above procedure. It only affects the numerical determination of M in a loop expansion. 



References 



[1] R. Efraty and V.P. Nair, Phys.Rev.Lett. 68 (1992) 2891; Phys.Rev. D 47 (1993) 
5601. 

[2] R. Pisarski, Physica A 158(1989) 246; Phys.Rev.Lett. 63, (1989) 1129; E. Braaten 
and R. Pisarski, Phys.Rev. D 42 (1990) 2156; Nucl.Phys. B 337 (1990) 569; ibid. B 
339, (1990) 310; Phys.Rev. D 45 (1992) 1827. 

[3] J. Frenkel and J.C. Taylor, Nucl.Phys. B 334 (1990) 199; J.C. Taylor and S.M.H. 
Wong, Nucl.Phys. B 346 (1990) 115. 

[4] R. Jackiw and V.P. Nair, Phys.Rev. D 48 (1993) 4991. 

[5] J.P. Blaizot and E. lancu, Phys.Rev.Lett. 70 (1993) 3376; Nucl.Phys. B 417 (1994) 
608. 

[6] V.P. Nair, Preprint CCNY-HEP 4/94. 

[7] A. Billoire, G. Lazarides and Q. Shafi, Phys. Lett. 103 B (1981) 450; O.K. Kalash- 
nikov, JETP Lett. 39 (1984) 405; T.S. Biro and B. MuUer, Nucl.Phys. A 561 



8 



(1993) 477; O. Philipsen, in Proceedings of the NATO workshop, Sintra, Portugal, 
March 1994 (to be pubhshed); W. BuchmuUer and O. Phihpsen, Preprint DESY-94- 
202 (November 1994). 

[8] J.M. Cornwall, Phys.Rev. D 26 (1982) 1453; J.M. Cornwall and J. Papavassiliou, 
Phys.Rev. D 40 (1989) 3474. 



9 



